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ABSTRACT: Black hole formation/evaporation in two-dimensional dilaton gravity 
can be described, in the limit where the number N of matter fields becomes large, by a 
set of second-order partial differential equations. In this paper we solve these equations 
nimierically. It is shown that, contrary to some previous suggestions, black holes evaporate 
completely a finite time after formation. A boundary condition is required to evolve the 
system beyond the naked singularity at the evaporation endpoint. It is argued that this 
may be naturally chosen so as to restore the system to the vacuum. The analysis also 
applies to the low-energy scattering of S'-wave fermions by four-dimensional extremal, 
magnetic, dilatonic black holes. 
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In recent work [CGHS], referred to as CGHS, a set of second-order partial differential 
equations were found to describe both 1 + 1 dimensional black hole formation/evaporation 
(including backreaction) and scattering of S'-wave matter by 3 + 1 dimensional extremal 
black holes. The CGHS equations apply to theories of gravity augmented by a scalar dilaton 
and N matter fields, and are the first term in a 1/A^ expansion of the full quantum theory. 
They provide a concrete arena in which to address the puzzles of quantum mechanical 
black holes. The CGHS equations apparently are not analytically soluble. Attempts have 
been made to remedy this problem by modifying the gravity-dilaton couplings. In [ABC] 
modifications were found which led to exactly soluble (even for finite N) quantum theories. 
Unfortunately these theories are unphysical because the Hamiltonian is unbounded from 
below [gs]. This difficulty was circumvented in [RST], henceforth referred to as RST, by 
imposing boundary conditions at the core of the black hole. These boundary conditions 
ruin the exact solubility of the full quantum theory, but not of the large- AT semiclassical 
equations. 

A fairly complete picture of black hole evaporation has emerged in the RST model. 
Black holes evaporate a finite time after formation, and the system can then be restored 
to the vacuum via an endpoint boundary condition [RST,AS]. Information is destroyed at 
the spacelike black hole singularity. 

While the results of RST are very enlightening, it is still important to try to understand 
the original CGHS equations. Large solubility was obtained in RST at the price of 
imposing an unnatural global symmetry. One may therefore wonder if the RST results are 
generic. (Indeed, previous analyses have suggested that solutions of the CGHS equations 
behave quite differently from those of the RST equations.) In fact we shall see in this 
paper that most, although not all, of the features of the RST model appear in a more 
generic setting. 

A number of previous attempts have been made to analyse the CGHS equations, but a 
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definitive conclusion has been surprisingly elusive. It was originally conjectured in CGHS 
that gravitational collapse did not lead to singularities. This was subsequently shown to 
be false [BDDO,ORST] in that black hole-like singularities do arise, and static solutions 
corresponding to a blackhole in equilibrium with a heat bath were found in [husk,nmr]. 
It was also shown [ORST] that if a large (apparent) black hole is formed, its apparent 
horizon shrinks due to Hawking evaporation. In [husk] it was suggested that the apparent 
horizon and singularity approach one another, meeting at future timelike infinity. In [thb] 
numerical evidence was given that they meet at a finite value of retarded time at future 
null infinity. This cataclysmic phenomena was referred to as a "thunderbolt" . 

In this paper we shall demonstrate numerically that in fact none of these final outcomes 
are realized. Instead, we find that the apparent horizon and singularity meet in a finite 
time, in agreement with results of Lowe [LOWE]. A boundary condition is needed at 
the junction point to continue the evolution of the system. We argue that this boundary 
condition can be chosen so that the system eventually settles back to the vacuum. However, 
because a spacelike singularity is present at intermediate stages, information is lost in the 
process of large N black hole formation/evaporation. 

Two-dimensional dilaton gravity coupled to N matter fields /j is described by the 
action 



where r = r ± u. In this gauge the "linear dilaton" vacuum solution may be written 




d^aV^ e-2<^(i? + 4(V(/>)2-f-4A2)--5^(V/,)2 . (actn) 



i=l 



We shall employ conformal gauge in which 



9++ = 9 — = 0, 



(Idv) 
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where here and henceforth we choose units so that A = 1. Classical black holes are formed 
by arbitrary infalling (i.e., d-f = 0) matter from I~ {t~ = — oo). The most general such 
classical solution is given by 



fd+f 



(csln) 



or a coordinate transformation thereof. The Penrose diagram for this process is illustrated 
in Figure 1. 

At the quantum level black holes evaporate. In the large N limit, with Ne^^ held fixed, 
the quantum theory is described by the semiclassical CGHS equations. A convenient form 
of the p — (j) equations is 

8Pd+d-(l) = -P\Ad+(t)d^(t) + A^e^^), 

{rpeq) 

2Pd+d-p = e-^'^{Ad+(f)d-(f) + A^e^^), 

where 

M (PPP) 

P' ^ -4e-4'^ + -e-'^t 
6 

These two equations imply the constraint equations, which are given by 

1 ^ TV 

e-2^(4a+0a+p - 2d\ct>) + - ^ d+hd+h - -{d+pd+p - d\p) + t+ = 0, itpp) 

i=l 

together with a similar equation for the ( ) components. t+{t-) is an arbitrary function 

of t'^{t~) determined by the initial data. 

In interpretation of solutions of (rpeq) it is useful to recall [CGHS,BDDO,dxbh] that 
1 + 1 dimensional dilaton gravity may be derived from S-waye reduction of 3+1 dimensional 
gravity. In that context 47re~^'^ measures the area of the two spheres of constant radius 
[dxbh]. Correspondingly a future (past) apparent horizon is a line along which d+(f){d-(f)) 
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vanishes [ORST]. (This is distinct from the event horizon, which is a null line from behind 
which one cannot escape to I"*".) This apparent horizon provides a useful two-dimensional 
definition of the boundary of a black hole. 

In the classical theory, the kinetic terms in the action degenerate at e"^*^ = 0, which 
may be thought of as the origin since the two spheres have zero area. In the large N 
quantum theory this degeneration occurs at the finite value (f)cr of (f) such that e~^^'='~ = 
as can be seen from inspection of (rpeq). Thus the quantum mechanical origin is at (per, 
and we shall not try to evolve the system beyond this point. 

We first compute the effects of sending in a shock wave of infalling matter along r"*" = 
with total energy M. Initial value data for this problem may consist of the values of p and 
(f) along two orthogonal null lines. We specify vacuum data along t+ = 0: 

p(0,T-) = 0, 

{vdta) 

0(O,T-) = -T-. 

Along I~ , the quantum equations exponentially approach the classical equations, and the 
initial data should agree with the classical solution (csln) with ^d+fd+f — Md{T~^): 

r+ > 0, 

T~ — OO, 



^ - - In 



iasdt) 



Of course it is not possible numerically to set initial data at r = — oo. Instead we specify 
the data at a large negative value Tq of t~; 

T+ > 0, 

00 = In M(l - e^^) + e^^-^o"] , ((^^^^) 



which induces an error of order e^o . We then check numerically for convergence at large 
negative Tq .'\ 

Our code approximates the spacetime using a null grid in r^. It computes the variables 
in some preset range < r"*" < t^^^ and then steps up in r~. The numerical computation 
cannot continue past a point where (f) = (l)cr, so it is instructed to revert to r"*" = and 
integrate along the next t~ step until the singularity is again encountered. In this way it 
fills in the region under (but not in the causal shadow of) the singularity. 

The results are plotted for M = .5 and A?" = 12 in figures 2a, 2b and 2c. Since N 
may be scaled out of the equations by shifting (f) the results are qualitatively independent 
of N. A black hole {i.e., a region where 5+0 > 0) forms a finite distance along the infall 
line. A spacelike singularity is produced when the Shockwave reaches the origin where 
(f) = (f)^j. {(f)^j. = for N = 12). All future- directed ingoing null geodesies inside the 
black hole terminate at the singularity. The apparent horizon becomes timelike due to 
Hawking radiation. Our computation clearly indicates that (contrary to previous claims) 
the apparent horizon meets the singularity in a finite proper time. At this "endpoint" 
denoted (r^, r^), the black hole has shrunk to zero size (as measured by e~^^ — ^). 

The endpoint singularity is visible to observers outside the black hole with coordinates 
> T^i so this constitutes a violation of cosmic censorship. Boundary conditions are 
required to evolve the system beyond the endpoint. Physically one expects that it should 
be possible to impose boundary conditions which restore the system to near-vacuum just 
after the endpoint. Indeed, inspection of Figure 3 shows that along the edge of the causal 
shadow of the singularity {r^ > t^jT^), the fields are very near their vacuum values. 
However, they are not exactly equal to their vacuum values. This is also expected: the 
outgoing Hawking radiation itself induces particle production. This leads to some flux 
crossing the line (r"*" > t^,t^) along which the black hole is seen to disappear, 
t Erroneous conclusions were reached in [thb] as a result of choosing a small value of Tq . 
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Potentially this infalling energy flux could recoUapse to form a new black hole. How- 
ever, we do not think that this is the case, because dilaton black holes (of all sizes) radiate 
energy at a rate of order N/A8'\. Therefore there is a threshold for black hole production: 
energy must be thrown in faster than it is radiated away. Such a threshold was found ana- 
lytically in the RST model, and numerical evidence for the existence of a similar threshold 
in the CGHS model is illustrated in figures 4 and 5. 

Although presumably no black hole is formed below the threshold, {i.e., 9+0 < 
everywhere), the CGHS equations are still ill-defined at </> = (per- However since this is 
a timelike line, it is presumably possible to choose reflecting boundary conditions so that 
the below-threshold radiation passes through the origin without collapsing. 

A plausible picture of black hole formation in the CGHS model is thus the following. 
The black hole evaporates completely in a flnite time and, with the help of an endpoint 
boundary condition, disappears completely. Weak, backscattered, infalling radiation re- 
mains, but this can be reflected through the origin. The excitations eventually dissipate 
and the system settles back down to the linear dilation vacuum. The Penrose diagram is 
illustrated in Figure 6. 

The CGHS model differs from the RST model in this last respect: In the case of RST, 
there is no backscattered radiation, and the system is identically vacuum in the causal 
future of the black hole endpoint. This appears to be a non- generic feature of the RST 
model, perhaps associated with the unnatural symmetries. 

We have also considered the semiclassical one-loop ghost-corrected equations derived 
in [fpbh], and given by (rpeq) with the flrst equation in (ppp) replaced by 



t The black hole mass M is also of order N with our conventions. 
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The extra term in (pfb) moves the critical values of into the complex plane and eliminates 
the singularity. The behavior of the system was analyzed in [fpbh]. In particular it was 
shown that the static solutions describe black holes with non-singular, asymptotically 
deSitter interiors in thermal equilibrium with a bath of Hawking radiation. A numerical 
simulation of shock wave collapse described by these equations is illustrated in Figure 7. A 
more detailed analysis is required to fully understand this system, but it appears that the 
shock wave creates a black hole which subsequently evaporates down to a stable remnant 
with a deSitter-like interior. It is apparently possible to create a remnant with an arbitrarily 
weak shock wave, which indicates that the remnant mass may be less than the linear dilaton 
vacuum. 
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FIGURE CAPTIONS 



Figure 1. Penrose diagram for classical black hole formation by an / shock wave with 
stress tensor ^d+fd+f = MS{t'^ — Tq). 

Figures 2a, 2b, 2c. Numerical simulation of black hole formation and evaporation. Initial 
conditions are specified along the left and lower boundaries of the plot corresponding to a 
null M — .5 shock wave along the left boundary. The contours in Figure 2a depict lines of 
constant (f) (the rippled dashes are an artifact of the plotting routine). The interior of the 
black hole is the region where these lines slope downward to the right, and the apparent 
horizon is the boundary of this region. The results depicted here depend on two numerical 
parameters: Tq , the initial t~ value, and dr, the numerical time step. In order to ensure 
that the results are not due to numerical errors we checked the convergence of our results 
as — > — oo and dr ^ 0. In Figure 2b we depict the singularity line (f) = (f)cr and the 
apparent horizon line d+cf) = for dr values ranging between 4 • 10"^ and 6.25 • IQ-^. It 
is evident that the curves converge. More quantitative tests show that the errors decrease 
like (ir^ and e^o , as expected. 

Figure 2b suggests that singularity and the horizon intersect. However, technically, our 
scheme cannot reach the intersection point (t^iT^) itself. The numerical scheme breaks 
down and it is impossible to continue the computation in the r"*" direction when cp = (per- 
Regardless of how fine our resolution is we can never reach the endpoint {t^,t^). There 
will always be a gap of one computation step dr between the last t~ value for which we 
calculated the horizon's position and the first that is singular. However, Figure 2c shows 
an enlargement of the region near (t^,t^). The different points correspond to the two 
finest (smallest dr) computations and the curves describe a Pade extrapolation through 
these points. The agreement between the computed points demonstrates the numerical 
convergence of our computation. The fact that the extrapolated curves intersect, at a 
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point just beyond the end of the computation, is clear evidence that {t^iT^) exist and 
are near (.3154, -.1093). 

Figure 3. Numerical simulation of black hole formation and evaporation. Initial condi- 
tions are specified along the left and lower boundaries of the plot corresponding to a null 
M = .5 shock wave along the left boundary. This plot extends to larger values of r"*" than 
those appearing in Figure 2a. It is evident that lines of constant cross into the causal 
shadow of the black hole endpoint along ~ —0.1. The nearly straight nature of those 
lines for large t'^ also indicates that the system is near its vacuum configuration far from 
the black hole. 

Figure 4. Energy is thrown towards the origin at a constant, above-threshhold rate for 
r"*" > 0, corresponding to a matter stress tensor ^d-\-fd-\-f — 0{t~^). A black hole forms 
and grows because matter accretion outpaces Hawking radiation. A spacelike singularity 
eventually crosses the null line t~ = where the numerical simulation ends. 

Figure 5. Energy is thrown towards the origin at a constant, below-threshhold rate for 
r"*" > 0, corresponding to a matter stress tensor ^d-\-fdj^f = .59{t'^). A black hole does 
not form because Hawking radiation outpaces matter accretion. No singularity crosses the 
null line t~ = where the numerical simulation ends. 

Figure 6. Conjectured Penrose diagram for large-N black hole formation/evaporation. 
The black hole disappears in a finite time. Boundary conditions can be chosen along the 
timelike segment of the line (f) = 4)cr after the endpoint in such a way that the system 
eventually settles back down to the vacuum. 

Figure 7. Numerical simulation of black hole formation and evaporation, as described 
by the one-loop ghost-corrected semiclassical equations with N = 12. Initial conditions 
are specified along the left and lower boundaries of the plot corresponding to a null M = 
.1 shock wave along the left boundary. The contours depict lines of constant 0. The 
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simulation is terminated when e^P diverges. This is apparently a coordinate divergence an 
infinite distance away from finite points in the spacetime, indicating the possible formation 
of a deSitter-like region inside the horizon. For large r"*" the system appears to settle down 
to a static remnant. 
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